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A fraction of AGN producing VHE 7-rays are located in galaxy clusters. The magnetic field 
present in the intra-cluster medium would lead to conversions of VlfE photons into axion-like par- 
ticles (ALPs), which are a generic prediction of several extensions of the Standard Model. ALPs 
produced in this way would traverse cosmological distances unaffected by the extragalactic back- 
ground light at variance with VHE photons which undergo a substantial absorption. Eventually, a 
nontrivial fraction of ALPs would re-convert into VHE photons in the magnetic field of the Milky 
Way. This mechanism produces a significant hardening of the VHE spectrum of AGN in galaxy 
clusters. As a specific example we consider the energy spectra of two observed VHE 7-ray sources 
located in galaxy clusters, namely lES 0414-1-009 at redshift 2 = 0.287 and Mkn 501 at z = 0.034. 
We find that the hardening in the observed spectra becomes relevant at _B > 1 TeV. The detec- 
tion of this signature would allow to indirectly probe the existence of ultra-light ALPs with mass 
rua ^ 10~* eV and photon- ALP coupling ga-y ^ 10"^" GeV~^ with the presently operating Imaging 
Atmospheric Cherenkov Telescopes like H.E.S.S., MAGIC, VERITAS and CANGAROO-III and 
even more likely with the planned detectors like CTA, HAWC and HiSCORE. An independent lab- 
oratory check of ultra-light ALPs invoked in this mechanism can be performed with the planned 
upgrade of the photon regeneration experiment ALPS at DESY and with the next generation solar 
axion detector lAXO. 

PACS numbers: 14.80. Mz , 95.85.Pw, PACS. 98.70Vc, 98.70.Rz, 98.35.Eg. 



I. INTRODUCTION 



Axion-like particles (ALPs) are very light pseudo-scalar bosons a with a two-photon coupling 077 which are pre- 
dicted by several extensions of the Standard Model like four-dimensional ordinary and supersymmetric models (see 
e.g. Jlljp); Kaluza-Klein theories (see e.g. Q) and especially superstring theories (see e.g. [3-01) (for a review, 
see |8|, |9|). In the presence of an external magnetic field, the 077 coupling entails that interaction eigenstates differ 
from propagation eigenstates thereby leading to the phenomenon of photon- ALP conversion 700 and in particular to 
photon- ALP oscillations [lol - [l2j . This mixing effect is exploited to search for generic ALPs in light-shining-through- 
the-wall experiments (see e.g. the ALPS 13] and the GammeV [iJl experiments), for solar ALPs (see e.g. the 
CAST experiment jtlSi]) and for ALP dark matter [13] in micro- wave cavity experiments (see e.g. the AD MX 
experiment [Tsj). 

Photon- ALP oscillations would also lead to intriguing signatures in astrophysical and cosmological observations [oj . 
In particular, over the last few years it has been realized that the 077 coupling can also produce detectable effects 
in the observations of distant active galactic nuclei (AGN), since photons emitted by these sources can mix with 
ALPs during their propagation through large-scale magnetic fields [l^. In this context, photon- ALP oscillations [l9l - 
[2H provide a natural mechanism to drastically reduce the absorption of very high-energy (VHE) photons by pair 
production processes (7vhe -I-7ebl ^ e+e") on the extragalactic background light (EBL) above roughly 100 GeV. In 
this respect, recent observations of cosmologically distant gamma-ray sources by ground-based gamma-ray Imaging 
Atmospheric Cherenkov Telescopes (lACTs) have revealed a surprising degree of transparency of the universe to 
VHE photons [13, [2^. This result has been confirmed by a recent systematic analysis performed on 25 VHE photon 
sources, which pointed out a suppression of the pair production during the propagation of VHE photons, named "pair 
production anomaly" [2i |. 

Oscillations between VHE photons and ALPs can represent an intriguing - but not unique, see e.g. (25l - [30| - 
possibility to explain this anomaly. Indeed, if VHE photons are transformed into a mixed photon-ALP state, the 
ALP component docs not suffer from absorption effects while it propagates and can therefore reach the Earth from 
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TABLE I: A sample of Blazars in clusters of galaxies [il. |45|]. 



Blazar 


Celestial coordinates 


redshift 


Mkn 501 


161i53m52.2s + 39d45m37s 


z = 0.034 


PKS 0548-322 


051i50m38.4s - 32dl6ml2.9s 


z = 0.069 


PKS 2005-489 


201i09m27.0s - 48d49m52s 


z = 0.071 


PKS 2155-304 


211i58m52.7s - 30dl3ml8s 


z = 0.116 


lES 1101-232 


Illi03m36.5s - 23d29m31s 


z = 0.186 


lES 0414+009 


041il6m52.96s + 01d05m24s 


z = 0.287 



distant sources even at very-high energies. In this sense, two complementary mechanisms have been proposed: a) 
VHE photon- ALP conversions in the magnetic fields around gamma-ray sources and then further back-conversions in 
the magnetic field of the Milky- Way [3l|, Is^l ; b) oscillations of VHE photons into ALPs in the turbulent extragalactic 
magnetic fields [ssl - lssj . Both these mechanisms are intriguing, but are affected by possible drawbacks. In particular, 
concerning the mechanism a) it is not clear at all whether a conversion 7 — > a actually takes place in all AGN mainly 
because their magnetic field is quite complicated and poorly known (see, e.g., (STj). On the other hand, mechanism 
b) requires intergalactic magnetic fields B 0.1 — 1 nG close to the current upper bounds, for which there is no firm 
observational evidence, or, alternatively, a very large 077 coupling, at a level already excluded by other independent 
observations.^ 

The aim of the present paper is to investigate in detail a third possibility concerning AGN hosted in clusters of 
galaxies, where a first conversion 7 — a occurs in the magnetic field of the cluster while a reconversion a ^ 7 happens 
in the magnetic field of the Milky Way. This scenario has two important advantages. First, we are dealing with 
magnetic fields that are known to a considerable extent. Second, since a sizable fraction of the original photons travel 
all the time as ALPs, the overall EBL absorption is strongly reduced. As a consequence, a significant hardening of 
the VHE spectra of distant AGN located in galaxy clusters is expected. 

The structure of the paper is as follows. In Sec. 2 we consider the observational evidence of AGN producing VHE 
photons in galaxy clusters. In Sec. 3 we review the formalism describing the photon- ALP mixing. In particular, we 
discuss the mechanism of photon- ALP oscillations in random magnetic fields which is relevant for 7 — a conversions 
in galaxy clusters. In Sec. 4 we address the effect of photon- ALP oscillations on VHE photon spectra, describing the 
7 — >■ a conversions in the intra-cluster magnetic fields, the absorption of photons in the intergalactic medium and 
the a — > 7 conversions in the Milky Way magnetic field. We show how our mechanism can lead to a characteristic 
hardening of VHE photon spectra for distant sources at i? > 10 TeV. As a specific example we discuss in Sec. 5 the 
effect on the energy spectra of two AGNs located in galaxy clusters, namely Mkn 501 at z = 0.034 |4(|| and lES 
0414-f009 at redshift z — 0.287 41]. Finally, in Sec. 6 we discuss future perspectives and we present our conclusions. 

II. ACTIVE GALACTIC NUCLEI IN GALAXY CLUSTERS 

The bulk of extra-galactic sources of VHE 7-rays has been identified with X-ray emitting AGN (see e.g. [l^]), which 
can be roughly described as an accretion disk around a supermassive black hole and two jets emanating from the 
centre and perpendicular to the disk. When one of the jet is directed towards us the AGN is called a Blazar. 

Blazars [4j] are divided into two broad groups: BL Lacs and Flat Spectrum Radio Quasars (FSRQs). BL Lacs are 
defined by the weakness of their thermal features like broad emission lines in their optical spectra. Accordingly, the 
nuclear region of BL Lacs, where the jet forms and accelerates is believed to be rather poor of soft photons. On the 
other hand, FSRQs display luminous broad emission lines, indicating the existence of photo-ionized clouds rapidly 
rotating around the central black hole and forming the so-called broad line region (BLR) at about one parsec from 
the centre. 

BL Lacs are believed to be FanarofF- Riley I type galaxies (46l - l48| . High spatial resolution optical observations 
indicate that the host galaxies of X-ray selected BL Lacs are typically giant elliptical galaxies [4§-[5l|. Imaging 
observations of the environment of 45 Blazars up to z = 0.65 indicate that the bulk of investigated objects avoid rich 



Recently, theoretical arguments have been discussed that the intergalactic medium is efficiently heated through generation of plasma- 
instabilities by powerful blazars [39( |. If this heating mechanism is at work, it would imply a model-dependent upper limit on the 
intergalactic magnetic field strength of _B ~ 10""^ nG. 
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clusters [52] ■ However, according to a recent study of Sloan Digitial Sky Survey data including spectroscopic redshift 
information for objects with z < 0.4, BL Lacs are found both in low and high density large scale environment [ssj . 
The sub-sample of VHE-emitting Blazars appears to populate both rich clusters as well as unspecific environments 
[S^ l- At this regard, in Table U we present a sample of Blazars at different redshifts which have been found to be 
located in galaxy clusters. Moreover, for a few Blazars not directly associated with a galaxy cluster, it can well happen 
that the line-of-sight traverses an intervening galaxy cluster. 

The existence of magnetic fields in galaxy clusters is well established through the observation of radio synchrotron 
emission as well as through the rotation measure of polarized radio sources (for a review, see e.g. [55|). The structure 
of the magnetic field within galaxy clusters has been subject to a number of studies, most notably Faraday rotation 
measurements and subsequent modeling of the underlying magnetic field structure [H, [s^. For well-observed objects 
like e.g. the Coma and Hydra A clusters, a Kolmogorov-type power spectrum has been found to fit the data for scales 
up to tens of kpc [13, with field strengths of a few up to ten fiG in the inner range (< 100 kpc) of the cluster. 
The cluster magnetic field connects very likely smoothly to magnetic fields on larger scales (filament), which may be 
close to nG strength according to simulations [H, 1131 . 

This motivates our choice of parameters throughout this paper: A field strength B = 1 fiG and a coherence length 
Ic = 10 kpc. According to the standard lore, the intra-cluster magnetic field B is modeled as a network of magnetic 
domains with a size equal to the coherence length. In every domain B is assumed to have the same strength but its 
direction is allowed to change randomly from one domain to another. The typical electron density in the intra-cluster 
medium is rie — 1.0 x 10~ cm" 



III. OSCILLATIONS OF PHOTONS INTO AXION-LIKE PARTICLES 




A. General considerations 

Photon- ALP mixing occurs in the presence of an external magnetic field B due to the interaction term 

Caj ^ -^ga^F^^F'"'a = ga^E-Ba , (1) 

where ga-y is the photon- ALP coupling constant (which has the dimension of an inverse energy). 

We consider throughout a monochromatic photon/ ALP beam of energy E propagating along the xs direction in a 
cold ionized and magnetized medium. It has been shown that for very relativistic ALPs and polarized photons, the 
beam propagation equation can be written in a Schrodinger-like form in which 2:3 takes the role of time (ll| 

- , (2) 

where Ai{x3) and ^2(2^3) are the photon linear polarization amplitudes along the xi and X2 axis, respectively, 0(^3) 
denotes the ALP amplitude and A4 represents the photon- ALP mixing matrix. We denote by T{x3, 0; E) the transfer 
function, namely the solution of Eq. ([2]) with initial condition T(0, Q;E) — 1. 

The mixing matrix A4 simplifies if we restrict our attention to the case in which B is homogeneous. We denote by 
By the transverse magnetic field, namely its component in the plane normal to the beam direction and we choose the 
X2-axis along By so that Bi vanishes. The linear photon polarization state parallel to the transverse field direction Bj 
is then denoted by and the orthogonal one by A±. Correspondingly, the mixing matrix can be written as (6^ l63j 

/ \ 

(3) 

^Bt/^ and Aa = —rn^/2E, where ma is the 
ALP mass. The term Api = —UJ^J2E accounts for plasma effects, where cjpi is the plasma frequency expressed as 
a function of the electron density in the medium rig as ojpi ~ 3.69 x 10""'^^ Y^nT/cnr^eV. The terms Ajy''^ describe 
the Cotton-Mouton effect, i.e. the birefringence of fluids in the presence of a transverse magnetic field. A vacuum 
Cotton-Mouton effect is expected from QED one-loop corrections to the photon polarization in the presence of an 
external magnetic field Aqed = |A™ — AjP^I cx _B|,, but this effect is completely negligible in the present context. 
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An off-diagonal Aji would induce the Faraday rotation, which is however totally irrelevant at VHE, and so it has 
been dropped. For the relevant parameters, we numerically find 

A,.-l.lxlO-(JL)"(^^)kpc-, 
Aqeo =.1-1x10- (4) (^)%pc-'. (4) 

For the above estimates that we will use in the following as benchmark values, we refer to the following physical 
inputs: The strength of _B-fields and the electron density Ue are typical values for galaxy clusters mentioned in 
Sectionini the value of the photon- ALP coupling is below the direct experimental bound g^^ < 8.8 x 10~^^ GeV~^ 
obtained by the CAST experiment for ma ^ 0.02 eV [l^, slightly better than the long-standing globular-cluster 
limit [63|. We recah that for ultra-light ALPs (rria < 10"^° eV) a more stringent limit gaj ^ I x 10~" GeV-^ ^ 
or even ga^ ^ 3 x 10~^^ GeV~^ has been derived from the absence of 7-rays from SN 1987A even if with a large 
uncertainty. 



B. Single magnetic domain 



Considering the propagation of photons in a single magnetic domain with a uniform B-field with Bi = 0, the 
component A± decouples away, and the propagation equations reduce to a 2-dimensional problem. Its solution 
follows from the diagonalization of the 2-dimensional mixing sub-matrix of A^o through a similarity transformation 
performed with an orthogonal matrix, parametrized by the rotation angle 9 which takes the value [111] 

9=larct'.n(-^^) . (5) 
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In particular, the probability for a photon emitted in the state to oscillate into an ALP after traveling a distance 
d is given by [ll| 



= sin220sin2 



^osc d> 



, sin^(Ao,cd/2 ) 
(Aoscd/2)2 



i^a,d r ;: -Z:' . (6) 



where the oscillation wave number is 

Aosc= [(A„-Apif -f4Ay'/' . (7) 
It proves useful to define a critical energy [l^ 

^ _ i?|A,-Api| 



9 A 



50|m2-w2| /io-6g\ /5x 10-"GeV-i 



(10-8eV)2 \ Bt J \ ga-y 
in terms of which the oscillation wave number can be rewritten as 



GeV 



Aosc = 2A,^^/1+ ( ^ ) . (8) 

From Eqs. ([S]) - (|5]) it follows that in the energy range E ^ Ec the photon- ALP mixing is maximal {9 ~ 7r/4) and the 
conversion probability becomes energy-independent. This is the so-called strong-mixing regime. Outside this regime 
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the conversion probability turns out to be energy- dependent and vanishingly small, so that Ec acquires the meaning 
of a low-energy cut-ofF. 

So far, we have been dealing with a beam containing polarized photons, but since at VHE the polarization cannot 
be measured we better assume that the beam is unpolarized. This is properly done by means of the polarization 
density matrix 



V a(x3) J 

which obeys the Liouville-Von Neumann equation ^] 

dx3 

associated with Eq. ([2])- Then it follows that the solution of Eq. pO)) is given by 

p(x3, E) = T(x3, 0; E) p(0) T^ix3,0; E) , 
where p(0) is the initial beam state. Note that for a uniform B even if we clearly have 



(9) 



T{x3„Q-E) = e'(E+Mo)x., 



Eq. (HH) reads 



p{x3,E) 



p(0)e 



-iMo X3 



(10) 

(11) 

(12) 
(13) 



C. Network of domains with random magnetic fields 



Since in the following we will consider VHE photons emitted by an AGN in a galaxy cluster, we have to deal with a 
more general situation than the one depicted in the previous Section. Indeed, as discussed in Sec. [Ill the intra-cluster 
B-field has a domain-like structure with size set by its coherence length. Although the strength of B is supposed to be 
the same in every domain its direction changes randomly from one domain to another. Therefore the propagation over 
many magnetic domains is clearly a truly 3-dimensional problem, because - due to the randomness of the direction 
of B - the mixing matrix A4 entering the beam propagation equation cannot be reduced to a block-diagonal form 
similar to Aia in all domains. Rather, we take the xi, X2, X3 coordinate system as fixed once and for all, and - 
denoting by ■(/'fc the angle between Bt and the X2 axis in the generic fc-th domain (1 < k < n) - we treat every ipk 
as a random variable in the range < f/'fc < 27r. During their path with a total length L in the galaxy cluster, the 
beam crosses n = L/lc domains, where Ic is the size of each domain: The set {i?fc}i<fc<n represents a given random 
realization of the beam propagation corresponding to the set {V'fc}i<fc<n- Accordingly, in each domain the matrix M 
takes the form [631 



Mk = 



' A^y Aaj sin-f/ife 

cos Ipk 

yAa^sinV'fc Aa^ COSIpk Aa 



with 



Aa;a; = A II sin^ Vfc + A_L cos^ Vfc 



(14) 



(15) 



(All 



I sin V'fe cos V'fc 



(16) 



Ayy — A» cos^ ipk + A± sin^ ^pk 



(17) 



Working in terms of the Eq. (|10p . after the propagation over n magnetic domains the density matrix is given by 
repeated use of Eq. (fT3| with Mq — j> Mk, namely [67] 



Pn =T{lpn, ■ ■ ■ ,^1) PoT^l/jn, ■ ■ ■ ,1pl) , 



(18) 
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where we have set 

n 
k=l 

with 

Tk = e"^''^ , (20) 

which is the transfer function in the fc-th domain (as explained above, the E'-dependent factor drops out from Eq. 
(IT51) and so it has been neglected). 

Since we do not know the particular configuration crossed by the beam during its propagation, in order to get an 
idea of the effect induced by the 7 — > a conversions it is useful to perform an ensemble average over all the possible 
realizations encompassing the l,...n domains. Assuming that the conversion probability in each magnetic 

domain [Eq. ([S])] is small, the average photon flux after n domains reads 

q = Pll + P22 = P^h-A ' (21) 
where is the emitted photon flux. For an initial unpolarized photon state 

pn = ^diag(l,l,0) (22) 

we obtain [g^ 

where the ALP-photon conversion probability Pi-l^ is provided by Eq. ^ . 



IV. THE MECHANISM 



Equipped with the results of the previous Section, we are now ready to investigate the implications of the mechanism 
we are considering for the conversions 7 — > a — > 7 in the VHE range. More specifically, we first address the 7—^0 
conversions in a galaxy cluster (Sec. IIV Ap . then the propagation of photons in the intergalactic medium where pair- 
production effects are relevant (Sec. IIVSP and finally the back-conversions a — 7 inside the Milky Way (Sec. lIVC)) . 

For definiteness, we will assume in the following a photon-ALP coupling ga-y ~ 5 x 10~^^ GeV""'^ and an ALP mass 
nia = 10~^ eV. With these values, the critical energy in Eq. (|8]) above which the strong-mixing regime in a galaxy 
cluster takes place is Ec ~ 50 GeV. We remark that for our input parameters, we can neglect 70a conversions in the 
intergalactic medium. Indeed, assuming a strength of the intergalactic magnetic field i? ^ I nG as required for sizable 
photon-ALP conversions, the critical energy is Ec — 500 TeV, far above the reach of the experimental detectability. 
Conversions in the intergalactic i?-fields would otherwise be relevant for values of nia ^ lO^^^'eV. Finally, since our 
goal is to understand the relevance of galaxy clusters for 7 ^> a conversion, we assume that no such conversion occurs 
inside the Blazar. 



A. Conversions in intra-cluster magnetic fields 



As discussed in Sec. |TT1 we assume a cellular structure for the intra-cluster magnetic fields, with domain size 
Ic — 10 kpc. Moreover, for rie ~ 1.0 x 10~'^cm~'^ the plasma frequency turns out to be ujp\ ~ 1.2 x 10~^^eV. With 
these input values we are in the strong-mixing regime, where the photon-ALP conversion probability is energy- 
independent. From Eqs. ([?!]) - the average flux of photons and ALPs coming out from the galaxy cluster are 
given by 



= p: 



CL 



7— )-7 7 

1^1 1 



(24) 
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FIG. 1: Evolution of the photon flux I-y and ALP flux la (the emitted photon flux is normahzed to 1) for the Blazar lES 
0414+009 at redshift z — 0.287 for a representative energy E = 8 TeV. The left panel represents the evolution inside the galaxy 
cluster, the central one the evolution in the intergalactic medium (photon absorption and ALP free propagation, with redshift 
on the horizontal axis), and the right one the evolution within the Milky Way. 

The evolution of these fluxes inside the galaxy cluster - with normalized to 1 - are represented in the left panel 
of Fig. [U considering the case of Blazar lES 0414+009 for a representative energy E = 8 TeV. It turns out that 
on average the fraction of VHE photons converting into ALPs within the cluster is ~ 30 %. However, significant 
variations are possible along a given line of sight, owing to the random nature of the intra-cluster magnetic field. In 
particular, it is straightforward to show that for an initially unpolarized photon beam the ALP flux coming out of 
the cluster can vary between 0% and 50%. 

B. Absorption of VHE photons on extragalactic background light 

The ALPs produced in a galaxy cluster propagate over cosmological distance undisturbed until they reach the 
Galaxy. On the other hand, VHE photons undergo absorption due to the pair-production process off EBL low energy 
photons 'j/'^HE^EBL _^ e+e^. The energy range 100 GeV ^ E ^ 10 TeV relevant for presently operating lACTs 
absorption is dominated by the interactions with optical/infrared EBL photons. The absorption rate for such a 



8 



process as a function of the incident VHE photon energy E is |68l-l7l 



r^{E) = / de —2— / ^a,^{P) , (25) 



ml/E 



where e and n™'^ are the EBL photon energy and number density, respectively, ^ is the cosine of the angle between the 
VHE and the EBL photon momenta and the limits of integration in both integrals are determined by the kinematic 
threshold of the process. Moreover 



cr^-y(/3) = 1.25 X 10^25 (1-/32-) 



2/3(/32-2) + (3-/34)logl±| 



cm^ , (26) 



is the cross-section for the pair-production process [7l|, UM ^ function of the electron velocity in the center-of-mass 
frame /3 = [1 - 2m^/£;e(l - ClY'^. 

Several realistic models for the EBL are available in the literature which rely upon different strategies (see, e.g., [tsI - 
Remarkably, they are basically in agreement with each other. Among the possible choices, we will employ as 
benchmark EBL the Franceschini-Rodighiero-Vaccari (FRV) model ,75, 76, 76]. We will also compare our results with 
the minimal EBL Kneiske model [8l|, [s^l, which provides a strict lower-limit for the attenuation of VHE 7-rays. As 
we will see, our conclusions will be rather independent of the choice of the EBL model, as the dominant part of the 
7-ray flux at Earth is mostly given by back-conversions a — 7, as it is already clear from our Fig. [1] 

According to conventional physics, for a given Blazar at distance do (corresponding to a redshift zq), the photon 
spectrum observed at Earth is given by 

/f(i?)=exp(-r^)/^;(i?o), (27) 

where I^{Eq) is the emitted spectrum with initial photon energy Eq = E{1 + zq) and is the optical depth which 
accounts for the EBL absorption and reads 



I dxT^{E,x) 
Jo 

dz T^{E,z) 



Ho Jo (l + z)v/f^A + a„(l + ^)3 



(28) 



where Hq — 73 km Mpc^^ s~^ is the Hubble constant, fim — 0.24 [83] is the matter density and 51a = 1 ^ is the 
dark energy density (assuming a flat cosmology). 

Instead in the presence of 7 — >■ a conversions in galaxy clusters, only the photons coming out of the cluster will be 
EBL absorbed, as shown in the central panel of Fig. [TJ Therefore, the average fluxes reaching the edge of the Milky 
Way turn out to be 

I^-^ = exp(-r,)/CL^exp(-r,)pCV° , 

= e = (l-P,%)/0 . (29) 

Given the strong attenuation of the photon flux due to the EBL absorption, we get I^'^ ^ 



C. Back-conversions in the Milky Way 



Ultimately the photon/ALP beam crosses the Milky Way before being detected. Observations over the last three 
decades have led to a rather detailed picture of the Milky Way magnetic field. Perhaps, its most important feature 
is that it consists in two components, a regular and a turbulent one. The latter component can be described by a 
cellular structure, with strength B ~ 1 x 10~^ G and domain size Imw — 10~^ pc [sj]. It is straightforward to realize 
that in this case the oscillation length Zqsc — 27r/Aosc is much larger than the domain size ^mWi so that the a — >■ 7 
conversion is vanishingly small. Therefore, in the following we restrict our attention to the regular component which 
is relevant for a ^ 7 conversions. 

Measurements of the Faraday rotation based on pulsar observations have shown that this component is parallel to 
the Galactic plane, apart from a possible dipole component at the galactic center and a small vertical component. Its 
strength varies between B ~ 1.4 x 10^^ G in the Solar neighborhood and B ~ 4.4 x 10^^ G in the inner Norma arm 
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|85l l8q |. Moreover, the associated radial coherence length is Ir — lOkpc [8g, while the vertical scale height is of the 
order of a few kpc. Inside the Milky Way disk the electron density is ~ 1.1 x 10~^ cm"^ [s^l, resulting in a plasma 
frequency Wpi ~ 4.1 x 10^^^ eV. 

The regular component of the galactic magnetic field is typically split into a disk and a halo part. The disk 
component is most often modeled as a logarithmic spiral with either symmetric or antisymmetric behavior with respect 
to the galactic plane (see e.g. [H, and references therein). Additionally, depending on whether the direction of the 
field in two different arms is the same or opposite, the model is called axisymmetric (ASS, also dubbed disymmetric, 
DSS, in [1^) or bisymmetric (BSS), respectively. Its typical scale height is of order 1 kpc. The halo component 
is instead typically taken as a purely toroidal (i.e. azimuthal) field [o^l, with a slightly larger vertical scale height, 
possibly extending up to > 3 kpc. The halo field can be different above and below the galactic plane. 

Several models have been proposed in the literature, based on combined analyses of Faraday rotation measurements 
of extra-galactic sources and of polarized galactic diffuse radio emission. However, statistical fits of these models to data 
are generally poor [8l,|9l|. Recently, evidence for an extra, out-of-plane, "X-shaped" component has been suggested 
jo^ . together with an indication of large vertical scale height of the order of 5 kpc for the halo field component. While 
we take the Jansson and Farrar model "o^] as our benchmark, we also consider the best fit model derived by Pshirkov 
et al. [9^ (see their table 3) from a global analysis of Faraday rotation measurements of extragalactic sources, to show 
the systematics associated to the poor knowledge of galactic magnetic fields. We further check some of the models 
discussed in [89'] , which give results similar to the ones of [93| . 

In order to compute the a — > 7 conversions in the Galaxy, we integrate Eq. (jlOp along each given galactic line of 
sight inside it. An illustrative sky map of the line-of-sight dependent probability for an ALP at the edge of the Galaxy 
to convert into a photon at Earth is shown in Fig. [5] for our chosen reference Jansson and Farrar magnetic field model 
(upper panel) and for the Pshirkov et al. model [93] (lower panel). The effect of the different magnetic field models 
is striking. In the case of Jansson and Farrar model the probability of a — > 7 conversion is generally larger than for 
the Pshirkov et al., especially towards the galactic center and in the southern galactic halo, due to the presence of 
the X-shaped field and to the large vertical scale height of the halo field. This would favor the detection of hardened 
VHE spectra from extragalactic sources if the mechanism we are proposing is indeed at work. 

The position of the two Blazars under consideration is also marked in Fig. [51 Both sources are located in a region 
where the back-conversion probability can be of the order of a few percent up to ^ 30%. We remark that the two 
sources considered here are located in directions where different magnetic field models give comparable results, within 
a factor of a few, for the conversion probability. 

Moreover, as we already explained, we show in the right panel of Fig. [T] the evolution of the photon and ALP fluxes 
along the line-of-sight to the source lES 0414+009 in the Milky- Way. The observable average VHE photon flux at 
Earth is then given by 

jE _ pMWjMW ^ /-[^ _ pMW\ jMW 



7 7— ^7 7 V 7^7 

exp (-r,) P^J,P,% + (1 - ^,^7) (1 - P,%)]l",- (30) 

We remind that this expression is valid only to get the average photon flux. In general, in order to calculate the 
observable photon flux at Earth one should propagate the full evolution matrix inside the cluster, the intergalactic 
medium and the Milky Way. 

Finally, we stress that the stochastic nature of the 7 — ;> a conversions in the magnetic held of a galaxy cluster 
implies that the flnal conversion probability can be considerably different from the average since the true magnetic 
held configuration along the line of sight is unknown. This fact entails that the photon fiux observed at Earth should 
be better characterized in terms of probability distribution functions, obtained by considering 7^-0 conversions 
over different realizations of the intra-cluster magnetic field. An example of these distributions is shown in Fig. [3] 
for the source lES 0414-1-009 at different energies. These distributions have been found by simulating 10^ different 
realizations of the intra-cluster magnetic field. For simplicity we have normalized the distributions to their maxima. 
In particular, the distributions at low and high energies reflect the distribution of Ij and la, respectively, at the exit of 
the cluster. At low energies {E ^ 100 GeV), since the conversion probability Pa-s-7 in the Milky Way is small (indeed, 
we are close to the critial energy Ec), the dominant contribution comes from the emitted - almost unabsorbed - 
photon flux. Conversely, at high energies {E > 10 TeV) almost all the emitted photons are absorbed and so the flnal 
photon flux is practically due to the reconverted ALPs in the Milky Way. At intermediate energies the distributions 
are a sort of "combination" between the emitted photon and ALP distributions. From Fig. |3] it is also obvious that 
the distributions are highly skewed with a changing asymmetry from low to high energies. We also must keep in mind 
that since photons of different energies cross the same (unknown) configuration of magnetic field, the distributions 
at different energies are not independent. The study of the correlations of the distributions at different energies is 
however beyond the scope of this work. 
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FIG. 2: Skymap in galactic coordinates of the a ^ 7 conversion probability, starting from a pure ALPs beam at the outside 
boundary of the Galaxy, for the Jansson and Farrar magnetic field model derived in [Q^l (upper panel) and the one of Pshirkov 
et al. [931 (lower panel), energy E = 8 TeV, coupling Qa-y = 5 x 10~^^ GeV~^ and rua = 10~^ eV. We also show the sky 
position of the two Blazars lES 0414+009 and Mkn 501. 

V. RESULTS 

After having described the details of the proposed mechanism, we discuss the observational signatures concerning 
the energy spectra of VHE photon sources. As specific examples we consider again the energy spectra of 1ES0414+009 
at redshift z — 0.287 and Mkn 501 aX z = 0.034. In order to show the effect of 7 a ^ 7 conversions on VHE 
photons, we exhibit in Fig.|4]the observable photon flux 1^ (normalized to the emitted one) as a function of energy for 
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FIG. 3: Probability distribution functions (normalized to the maximum) for the observable photon flux I., from 1ES0414+009 
at different energies. 

lES 0414+009 and for Mkn 501 for both the models of Jansson and Farrar .92] (upper plot) and Pshirkov et al. (93j 
(lower plot). The solid black curves represent the flux expected in the presence of EBL absorption only. According 
to conventional physics, it turns out that the flux gets dramatically suppressed at high energies {E > 1 TeV), the 
farthest the source the lowest the energy. 

Including the effect of 7 — > a ^ 7 conversions we see that the photon flux at high energy gets strongly enhanced 
with respect to the expectation in the presence of conventional physics. In particular, the continuous red curves 
represent the photon flux in the presence of 7 a 7 conversions, averaged over many realizations of the intra- 
cluster magnetic field. Indeed, the effect is striking. Since the photon- ALP conversion probability in the strong-mixing 
regime in which we are working is energy-independent, the photon flux displays a plateau - instead of a sharp drop - 
at which its intensity depends on the adopted Galactic magnetic field model. On average it turns out to be between 
10 % and 3 % of the emitted value for lES 0414-^009, and between 10 % and 1 % for Mkn 501. Therefore, the existence 
of 7 — >■ a — > 7 conversions produces a considerable hardening of the spectrum at high enough energies, thereby making 
it possible to detect VHE photons in a range where no observable signal would be expected according to conventional 
physics. 

However, we remark that due to the stochastic nature of the 7 — ?> a conversions in the random intra-cluster magnetic 
fields, the observable photon flux could show large fluctuations depending on the realization of the magnetic network 
crossed during the propagation. In view of these fluctuations, the systematic effect of our poor knowledge of Galactic 
magnetic fields plays a relatively minor role (as demonstrated in Fig.S]), although detailed observations of hard spectra 
in the VHE regime could be sensitive to the structure of the Galactic magnetic field. 
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FIG. 4: Observable photon flux (the emitted photon flux is normalized to 1) as a function of energy for lES 0414+009 and 
for Mkn 501 using the Jansson and Farrar model ^93] (upper plot) and the Pshirkov et al. model ^] (lower plot) respectively. 
The black solid curve represents the flux expected in the presence of only absorption onto EBL. We are using the FRV EBL 
model [T^, [T^. The solid red curve represents the average photon flux in the presence of 7 — >■ a — >■ 7 conversions. The dashed 
curve corresponds to conversions for a particular realization of the intra-cluster magnetic field. The shaded band is the envelope 
of the results on all the possible realizations of the intra-cluster magnetic field. 
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TABLE II: Fit parameters for the spectra of lES 0414+009 and Mkn 501 corrected with the FRV EBL model. The spectra are 
described with a power law, diV/d_E = No{E/Eoy^ and a logarithmic parabola dN/dE = 7Vo(£'/£o)'^"'""'^''^°'- The de- 
correlation energies are Eo = 0.27 (0.28) TeV and Eo = 1.98 TeV for lES 0414+009 (with ALPs) and Mkn 501 (both scenarios), 
respectively, and are held fixed. The normalization A^o is given in units of TeV~'^cm~^s~^. 





lES 0414+009 


Mkn 


501 


Fit parameters 


w/o ALPs w ALPs 


w/o ALPs 


w ALPs 


Normalization No x 10~^^ 


2.03 + 0.24 1.71 + 0.21 


4.60 + 0.06 


3.25 + 0.04 


Power-law index F 


2.32 + 0.30 1.86 + 0.39 


2.11+0.02 


2.12 + 0.02 


Curvature /3 




0.16 + 0.02 


0.10 + 0.02 




0.23 0.60 


1.36 


1.94 


d.o.f. 


4 4 


14 


14 



An example of a particular realization is shown by the dashed curve. In this specific case we see that the observable 
photon flux at high energies can be even larger than the average one. However, if one considers many realizations of 
the intra-cluster magnetic fields one obtains the shaded band as envelope of the results. Therefore, depending on the 
particular magnetic realization crossed by the photons, it is also possible to have cases in which the suppression of 
the photon flux is stronger than in the presence of conventional physics. Nevertheless, at high energies from Fig. 3 
one infers that the cases in which 1^ is enhanced are more probable. In general, 7 — > a — > 7 conversions cannot be 
regarded as a universal mechanism to produce an hardening in the spectrum of VHE photons. Conversely, if a large 
number of AGN in clusters were observed with large spectral variations this would be a supporting indication for the 
presence of ALPs. 

In Fig. E] we consider the spectra of lES 0414+009 (upper plot) and Mkn 501 (lower plot) using the FRV EBL 
model 76]. The observed spectra (black bullets) are corrected for EBL absorption with and without ALPs. This 
yields the data points shown as dark and light gray bullets, respectively. The EBL absorption corrected spectra are 
subsequently fitted with a power law (lES 0414+009) and a logarithmic parabola (Mkn 501), represented by the light 
blue (no ALPs) and dark blue bow ties (with ALPs). The final fit parameters are shown in Table HIl The difference 
in the observed spectra becomes prominent beyond ~ 4 TeV (lES 0414+009) and ^ 30 TeV (Mkn 501), respectively, 
where the 7 a — >■ 7 conversion (dark blue dashed lines) predicts a substantial flux enhancement. Unfortunately, with 
the present data it is impossible to distinguish the two scenarios. We also remark that with ALPs, the emitted photon 
flux at low energies has to be larger than in the case of conventional physics in order to describe the observations. 
The reason is that more photons convert into ALPs in the cluster than the other way around in the Milky Way and 
the absorption at low energies is negligible. 

The same analysis is repeated for the minimal EBL Kneiske model [Sll, 113 and the corresponding results are shown 
in Fig. [SI The difference between the EBL models is marginal since the spectra are dominated at high energies by 
the ALPs that re-convert into photons in the Milky Way. 

Finally, the boost factor of the flux intensity with respect to conventional physics 



is plotted in the lower panels Fig. [5H6]for the FRV EBL model [75|,l26[ and for the minimal EBL Kneiske model |8ll 
respectively. One realizes that when ALP-photon conversions are effective, they can produce B ^ 10^ thereby implying 
a strong enhancement of the observable photon flux with respect to expectations based on conventional physics alone. 



VI. CONCLUSIONS 



VHE 7-ray observations offer the possibility to indirectly probe the existence of ALPs predicted in many extensions 
of the Standard Model. In this respect, 7 -> a — 7 conversions of VHE photons in the presence of cosmic magnetic 
fields have been recently proposed as an intriguing mechanism to explain the surprising high degree of transparency 
of the Universe to VHE photons recently observed in different high-redshift sources. In the present work we have 
explored the further possibility of 7 a conversions in the magnetic field of galaxy clusters which frequently host 
VHE emitting Blazars and subsequent a ^ j regeneration in the magnetic field of the Galaxy. We have shown that 
this mechanism can produce a significant hardening of the VHE photon spectrum of Blazars located in clusters of 
galaxies. More specifically, the signature of this effect in lACTs allows to infer the existence of ultra-light ALPs 
with mass TOq ^ 10~* eV and photon- ALP coupling g^-y ^ 10~^" GeV~^. We expect that such a signature would 
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FIG. 5: The spectra of lES 0414+009 t41|] (upper plot) and Mkn 501 [401 (lower plot). We are using the FRV EBL 
model [t^, [t^. The upper panels show the observed spectra corrected for absorption with and without ALPs along with the 
extrapolated spectra, whereas the lower panels depict axion boost factor B, i.e. the ratio of the observed fluxes with and 
without ALPs. 
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FIG. 6: Same as Fig. [5] but with the minimal EBL Kneiske model [8ll.l8 
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start to emerge at energies E ^ 1 TeV. Therefore, they are barely testable with the present generation of lACTs like 
H.E.S.S. [al], MAGIC (95i], VERITAS and CANGAROO-III covering energies in the range from ~ 50 GeV 
to ~ 50 TeV. A clear-cut check can only come from the planned Cherenkov Telescope Array (CTA) and High 
Altitude Water Cherenkov Experiment (HAWC) [11], rea chin g energies of 100 TeV with much higher sensitivity, or 
even with the HiSCORE detector, reaching PeV energies jlOO| . 

A peculiar feature of the proposed mechanism is that since the magnetic fields in galaxy clusters have a turbulent 
nature, the 7 — a conversion probabilities present a large variance depending on the line of sight crossed by the 
photon/ALP beam inside the cluster. This fact suggests that our proposal could not be an universal mechanism to 
produce the transparency of the Universe to VHE photons. But if a large number of AGN in galaxy clusters were 
observed with large spectral variations this would be a positive evidence for our proposal. 

We also remark that even if an AGN is not located inside a galaxy cluster, there is a nontrivial chance that in some 
cases its line of sight crosses a cluster of galaxies. We plan to investigate what happens in this instance in a future 
work. 

Remarkably, an independent laboratory check of the ultra-light ALPs discussed in our scenario can be performed 
with the planned upgrade of the photon regeneration experiment ALPS at DESY [l3l | and with the next generation 
solar axion detector lAXO (International Axion Observatory) TOlj . This confirms once more the nice synergy between 
astrophysical and laboratory searches to find axion-like particles. 
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